Implementation of multi-walker quantum walks with cavity grid 
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We show how multi-walker quantum walks can be implemented in a quantum quincunx created via 
cavity quantum electrodynamics. The implementation of a quantum walk with a multi-walker opens 
up the interesting possibility to introduce entanglement and more advanced walks. With different 
coin tosses and initial states the multi- walker quantum walk shows different probability distributions 
which deviate strongly from the classical random walks with quadratic enhanced spreadings and 
localization effects. By introducing decoherence, the transition from quantum walks to the classical 
versions is observed. We introduce the average fidelity decay as a signature to investigate the 
decoherence-induced irreversibility of quantum walks. 

PACS numbers: 03.67.Ac, 42.50.Pq, 74.50.+r 



I. INTRODUCTION 

Quantum walks (QWs) Qj offer an alternative ap- 
proach to implement quantum algorithms [2( compared 
to the typical circuit-based 0] or measurement based Q 
models for quantum algorithms. There are two types of 
QWs: the continuous time QWs [f| and the discrete time 
QWs [f| ■ The realization for QWs have been proposed in 
quantum optics [7} , ion t rap [|| |9| , cavity quantum elec- 
trodynamics (QED) [l0|, [H[ and optical lattice [12| sys- 
tems in a decade. QWs on single photons [13], trapped 
ions [14| and neutral atoms [151 ] have been realized in 
the laboratory, respectively. More recently, higher di- 
mensional QWs [16| and QWs involving more particles 
[17H24| are studied. This reveals the additional features 
offere d by quantum mechanics, such as quantum correla- 
tions |24| and indistinguishability Compared to the pro- 
posals on QWs with single- walker, we extend the QWs 
by using multi-walker. The implementation of a QW 
with a multi-walker opens up the interesting possibility 
to introduce entanglement and more advanced walks. 

To create a QW, the following steps are followed: a 
particular QW is chosen, in our case a discrete multi- 
walker QW with joint coins each of which decides to the 
corresponding walker's position shifts. Furthermore, a 
signature for QW behavior is identified, such as enhanced 
diffusion or uniformity of the distribution for the walk- 
ers' degree of freedom. Then a physical system is cho- 
sen whose Hamiltonian dynamics match the evolution of 
the QW. Finally, open system dynamics are incorporated 
into the analysis in order to account for non-unitary evo- 
lution as well as to incorporate realistic measurement into 
the model. 

Compared to random walks (RWs), QWs are re- 
versible. The irreversibility due to decoherence trans- 
mits the QW to RW . It is of great interest to show the 
variation of the irreversibility in the time evolution of the 
QW. The probability distribution and standard deviation 
of the distribution are used to study the irreversibility 
of QWs in the present of decoherence. However those 
approaches are not operational. That means neither 
of them provides a way for direct experimental obser- 



vation of the irreversibility of the QW. In the present 
of decoherence, except for the probability distribution 
and the standard deviation, we introduce an operational 
measurement — the average fidelity decay (AFD) j25[|26|. 
The AFD can reveal the response of the system to the 
decoherence which is closely related to the properties of 
both the system and the environment, and provide an 
experimentally available way to monitor the detrimental 
influence on the QW with different decoherence sources. 



II. QUANTUM WALKS WITH ONE- AND 
MULTI-WALKER 



Let us first briefly review a single-walker QW. The 
Hilbert space of the walker+coin is given by a tensor 
product 



H=H W <S)H C 



(1) 



of the walker space % w and the two-dimensional (2D) 
coin space 



Wc = span{|-l),|l)}. 



(2) 



We consider a walker starting the QW from the origin, 
i.e., the initial state has the form 



IV^ini = IV'O), 



(3) 



where |V>o) w an( i IV') C denote the initial state of the walker 
and coin respectively. After TV steps of the QW, the state 
of the walker+coin is given by 



hM = u N |^) ini 



(4) 



= 5>_i(j, N) |Vi) w 1-1) +Pi(j, N) |^) w |1> , 
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where the unitary propagator U has the form 

U = S{l®C). (5) 



The probability distribution generated by the QW is 
given by 
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The coin operator C flips the state of the coin before the 
walker is displaced. In principle, C can be an arbitrary 
unitary operation on the coin space Ji c . We choose the 
most studied case of the Hadamard coin, denoted by H = 

/y/2, which is defined by its action on the basis 



1 -I 

states. 



H\±1) = (\-1)T\1) 



(7) 



After the coin flip, the step operator S displaces the 
walker from its current state according to its coin state 



shk>j±i>— H^±*>J±i> 



(8) 



where 6 is the step size. The coefficients p±i(J, N) repre- 
sent the probability amplitudes of finding the walker at 
\tpj) after N steps of the QW with the coin state |±1). 

As an extension, we choose the ./V-walker QW over 
circles in phase space, which arises naturally for TV har- 
monic oscillators. Points in phase space correspond to 
the oscillator position-momentum pair (x,p), which we 
henceforth refer to as the phase space 'location'. 

For the discrete N- walker QW on the circles, each of 
the walker's location as a point in phase space is replaced 
by a localized wavefunction centered at location (x,p), 
and the random flips are replaced by joint quantum coins 
given by qubits, which are flipped by a unitary operation 
and then entangled with the oscillators by free evolution. 
An example of TV-walker state is the product state 



'■>N 



(9) 



An example of the coin state is | ip) c = | c\ , c 2 , . . . , cn ) ■ The 
archetypal discrete time iV-walker QW consists of two 
building blocks — a coin operator C and a step operator 
S. The coin is essentially an ancillary parameter that is 
used by the step operator to decide how to propagate the 
walker. The simplest example of the coin operation is to 
apply a Hadamard transformation to each qubit in the 
decomposition of Eq. ([9]). This internal transformation 
is separable, in the sense that it does not produce entan- 
glement between the spatial degree of freedom. Other 
choices for the coin operations include the entangling coin 
operation, the discrete Fourier transform (DFT) and the 
Grover operation. A common choice of step operator is 



= |0i + c x 8) ® |02 + c 2 S) 
where c, G { — 1, 1} and 

I0i> : 



(10) 
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(ii) 



for j = 1, ..., N is the phase state. 

Using two-walker QW as an example, one obvious 
generalization of coin tossing operations is to apply a 
Hadamard transformation H to each qubit of the coin 
state. This choice can be viewed as two independent 
coin tosses on each qubit of the coin state. The transfor- 
mation is 



Ci=H®H 
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(12) 



This coin operation is separable, in the sense that it does 
not produce entanglement between the spatial degrees of 
freedom. 

In principle, any unitary transformations on the coin 
state can replace Hadamard transformation and be used 
as coin tosses. Now we introduce non-separated coin 
tosses. One obvious generalization of coin tosses, which 
is not separable and does produce entanglement between 
the coin qubits, is the root of SWAP gate operation, de- 
fined as follows 



C 2 = ViSWAP = 



a o o o\ 
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(13) 



Another choice of non-separated ncoin tosses — the dis- 
crete Fourier transform (DFT) D, defined as follows: 



C» = D = - 




(14) 



transforms any coin translation eigenstate into an equally 
weighted superposition of all the eigenstates and entan- 
gles the spatial degrees of freedom of coin qubits. 

There arc, of course, an infinite variety of other non- 
scparablc choices for the coin tosses by employing differ- 
ent phase relationships. Finally we introduce the Grover 
operator as a non-separated coin toss defined as follows 
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(15) 



We now investigate the effect of measurement on a 
two-walker system. After fcth steps of two-walker QW, 
the state of the walker+coin system becomes \ipk) = 
(SC) k \ipini). Suppose we perform a measurement and 
detect the first (second) walker at position <fi. Then the 
state will be projected into 



-Pi(2)(</0 = (0|Tr 1(2) p w 



(16) 



where p w = Tr c (\ip k ) (4>k\)- 
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FIG. 1: (Color online.) The ln-ln plot of the position spreads 
<t(0) of the first walker taking the ideal two- walkers QWs with 
different coin tosses (a) DFT D, (b) H®H, (c) ViSWAP and 
(d) Grover G as functions of the position cj> for different initial 
coin states |*) cl = (|1) + |-1) )® ( |1> + |-1> )/2 (♦), |*) c2 = 
|1)®|1) (■) and |*> c3 = (|l) + i|-l>)®(|l>+i|-l))/2 
(*)• 



s± As 


Wcl 


Wc2 


McS 


D 


0.940 ± 0.003 


0.880 ± 0.002 


0.887 ± 0.003 


H®H 


0.979 ± 0.005 


0.897 ± 0.007 


0.994 ± 0.004 


ViSWAP 


0.965 ± 0.001 





0.965 ± 0.001 


G 





0.911 ±0.002 


0.973 ± 0.006 



TABLE I: Linear regression results for ln<r(0) = (s±As) In A^ 
for the ideal case with the step size S = 0.8 in ln-ln scale shown 
in Fig. 1. 



the coin. The time dependence of the standard deviation 
for QW is plotted in Fig. 1 and the corresponding slopes 
Aa/ At (In a (<j))/ InN) are presented in Table I. We ob- 
serve that the standard deviation for QW with the same 
coin toss depends on the symmetry of the initial states. 
However the quantum behavior can always be observed 
except for the two cases — the -y/iSWAP QW with the 
initial coin state |1) <3 |1) and the Grover G QW with 
the initial coin state (|1) + |-1))®(|1) + |-1)) /2, which 
shows the localization effect on QW where the walker's 
spread becomes constant [28|, [29| ■ 

Any real implementation of a quantum system must 
deal with the issue of decoherence, which tends to de- 
stroy quantum correlations. The entanglement decay due 
to noise leads the transition from QWs to RWs. In the 
present of decoherence, except for the probability dis- 
tribution and the standard deviation of the distribution 
another operational measurement here is introduced as a 
new signature for QWs. We measure the irreversibility of 
the QW due to decoherence by the AFD which is defined 
based on fidelity decay, that is the square modulus of the 
overlap between two time dependent final states with the 
same initial state under the time evolution without and 
with decoherence respectively. The AFD [26[ is defined 



AFD = Tr[p(t)\V N )(* N \] t 



(18) 



where p(t) is the density matrix of the walker+coin sys- 
tem after the time evolution. The AFD reveals the re- 
sponse of the walker+coin system to the decoherence and 
is sensitive to different sources of decoherence. 



The dispersion of the distribution (1161) is especially 
important. As moments are not particulary useful for 
distributions over compact domains, other strategies are 
needed. For the phase distribution over the domain 
[0, 27r], Holevo's version [27[ of the standard deviation 



a{4>) 



d(j>P{4>)e^ 



(17) 



is particulary useful as it reduces to the ordinary stan- 
dard deviation for small spreads and is sensible when the 
dispersion is large over the domain. 

We begin our analysis the results of two-walker QW 
with different coin tosses. The initial conditions for the 
coin state were chosen to be the separable state com- 
posed of all qubits in the states |1), ( |1) + | — 1) )/\/2, and 
( |1) + i | — 1) ) /V%, which lead to three different proba- 
bility distributions. 

For the case of separable transformation with separable 
initial conditions, the different walkers behave indepen- 
dently; thus, the variance can be expressed in terms of 
one- walker case. Furthermore, a bias could be introduced 
into the transformation and give a different weighted su- 
perposition of translation cigenstatcs. Thus the phase 
distribution of walker also depends on the initial state of 



III. IMPLEMENTATION OF MULTI- WALKER 
QWS ON CIRCLES 

Multi-walker QWs can be implemented with scalable 
cavity grid [3(| with superconducting circuits [U |32 |. 
The walkers arc represented by the cavity modes and the 
coin states {| — 1) , |1)} are encoded with charge qubits. 
The cavity grid consists of cavity modes belonging to 
./Vh horizontal (H) and Ny vertical (V) cavities, H cav = 
Yj-j^iufajCLj + Yli=i u Ybjbj' coupled to one charge 



^=i w j "yj^ ^j=i w j "j"j> 
qubit at each intersection (i, j) 

chitecture: 



generalizing to a 2D ar- 



H cav . qb = Y. 1% < x l \a"^ + °<) + alih + b)) 



1,3 



H(V) 



(19) 

The coupling gf- y " ; between the horizontal (vertical) 
cavity mode i (j) and the charge qubits is switchable 
by the external electric field on each charge qubit J33| . 
Eq. (|19|) leads to the Jaynes-Cummings (JC) model and 
the cavity-mediated interaction between qubits. 

For a multi-walker QW, one can implement the coin 
operation on charge qubits capacitively coupled to a ver- 
tical cavity mode (j) with cavity-assisted interaction. 



The conditional position shifts of each walker can be im- 
plemented for each horizontal cavity mode i under the 
free evolution of the JC interactions between horizontal 
cavity mode i and charge qubit (i, j). 

The multi-walker QW with cavity grid can be imple- 
mented by steps as follows. Step I, the state of the cavity 
grid is prepared in a certain initial state. Step II, the cou- 
plings between the vertical cavity mode and the charge 
qubits arc turned on to implement the coin operation on 
charge qubits with cavity-assisted interactions. Step III, 
the couplings between the vertical cavity mode and the 
charge qubits are turned off and those between the hor- 
izontal cavity modes and the charge qubits are turned 
on to implement the conditional position shifts of each 
walker due to the coin states. Then, we repeat steps II 
and III for the next step of multi- walker QW. 



A. The Conditional Phase Shifts 

Using a two-walker QW as an example, we consider a 
system including two two-level charge qubit coupled to a 
cavity grid with the structure mentioned above. The cou- 
pling between each charge qubit and the corresponding 
horizontal cavity field is used to implement conditional 
phase shifts on circles and the charge states are used to 
implement quantum coins, each of them with two possible 
operations. The physical implementation of the condi- 
tional phase shifts of multi- walker quantum walks can be 
implemented by the free-evolution of the cavity-assisted 
interaction Eq. ([T5|. which can be rewritten in the JC 
model [U 

Hjc = E [wcataj + y u{ + g(a)al + a^)] (20) 

.7=1,2 

with w a and ui c the coin and cavity frequencies, respec- 
tively, and g the coupling strength. In the dispersive 
regime, 



|A| = |« a -w c | > .g, 



(21) 



and in a rotating frame, the effective interaction Hamil- 
tonian is 



.7=1,2 

with the cavity pull of the resonator 
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<£_ 

A 



(22) 



(23) 



After the coin flipping operation, the time evolutions of 
the interactive Hamiltonian give the conditional position 
shifts due to the charge states 



Uj = expiiAOjolajai), 



(24) 



for j = 1, 2, where A6ii 2 ) is the size of one step for walker 
1(2) which depends on the effective coupling \ and the 
time duration £1(2)- 



B. The Coin Tosses 

1. The Separable Coin Toss 

The coin tosses can be implemented on charge qubits 
with cavity-assisted interactions. Now we turn off the 
couplings between the charge qubits and the horizon- 
tal cavities. Then we turn on the coupling between the 
charge qubits and the vertical cavity subsequently. The 
time-dependent driving field applying on the vertical cav- 
ity [IH 



Ha 



#)( 



£t e -«"dt + y^t 



(25) 



can be used to implement the separable coin toss. It is 
sufficient to let e(t) be a square wave so e is a constant 
(e = when the field is off). In the dispersive regime 
and in a frame rotating at Wd for the qubit and the res- 
onator, Hj = iJj c + Hd can be replaced by the effective 
Hamiltonian 



S d 



O r 



H> = X tfbai - -£di + -^H - 6 dc tfb + e(&t + b) 



with 



Wd 



u d 



o - 2ge 

"R — ~i — 
Ode 



(26) 



(27) 



(28) 



the Rabi frequency. 

The first term in the above expression effects the coin- 
induced walker phase shift. The atom transition is an 
ac-Stark shifted by g 2 wb/S.. To implement 



H®H = n j= i i2 cxp [iiH^Roi/2] 



on the coin, we choose 
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(29) 



(30) 



with pulse duration tu = 7r/2fiR. 

After the coin flipping we shut off both the external 
field and the coupling between the charge qubits and the 
vertical cavity, and turn on the coupling between the 
charge qubits and the horizontal cavities. The free evo- 
lution continues for a duration tj for charge qubit j for 
the conditional phase shifts. 



2. The ViSWAP Coin Toss 

If one turn on the couplings between the two charge 
qubits and the vertical cavity at the same time, the ef- 
fective Hamiltonian of Hjc can be written as 

1 , 



H 
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X(ai.di + aia%' 



(31) 



The forth term is the induced dipole-dipolc interaction 
between the two charge qubits, which can be used to 
implement the \/— zSWAP coin toss on the charge states 
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2q = CXp 



*{M+b E^ 



ViSWAP. (33) 



3. The DFT Coin Toss 

The similar system can be used to implement a 2D 
discrete Fourier transform (DFT) coin toss, defined in 
Eq. (|14p . Note that the Hadamard transformation is the 
ID DFT [H. 

For the coin flipping, if we choose the duration tj to 
increase geometrically with qubit number as xtj = 2 J 7r/4, 
then the Hamiltonian Hf nt = xwba{ generate a unitary 
transformation 



D 



T^'bfs 



exp 



(34) 



where the electronic operator T provides a binary order- 
ing of the qubits: 



.7=1,2 



(35) 



The eigenvectors of the operator T are the electronic 
number states 



f = j>|fc) 



(36) 



fc=0 



where \k) = \S 



n, 



-tJV-1 



N 



IS 



andfc = S^ _1 + Sjv-ix2 



JV-l, 



JV-l 



N-2 



- O ISfi)!, Sj =0,1 

-S 1 x2°. The binary 
expansion for k is thus just the string SnSn-i-.-Si. 

This unitary operation can be used to implement a 2D 
DFT shown in Eq. (|14[) . First we turn on the coupling 
between charge qubit 1 and the vertical cavity and after 
time duration t\ it is turned off. Then we turn on that 
between charge qubit 2 and the vertical cavity and after 
time duration t? it is turned off. Thus the DFT operation 
on two charge qubits is realized. 



4- The Grover Coin Toss 

The induced dipole-dipolc interaction between two 
charge qubits can be used to implement the Grover coin 



toss [35J. Now we turn on the coupling between the 
charge qubits and the vertical cavity at the same time 
and apply a drive field (|2"5)) on the vertical cavity. The 
effective Hamiltonian of the system in the dispersive limit 
becomes 



^ q = E (* hrba i 



0/ 



i=i,2 



+ X(&+°1 + v-°l) - $dM + eft + b). (37) 

If JIr ^> <5da,<5dr,5, we can get the time evolution of the 
system in the interaction picture 



um 



-iH a t-iH c t 



(38) 



with 



H = R /2 Y, H ( 39 ) 

J=l,2 
He = X(0"l^2 + Ol^j) + h - C - 

For choosing \t = it/8 and SIr/x = 16m + 4 for m an 
integer, we can get 



Ui(t) = -G. 



(40) 



So by choosing appropriate values of parameters, we can 
generate a two-qubit Grover operation on two charge 
qubits. 

The multi-walker QWs can be implemented by three 
steps as follows. First, the cavity grid is prepared in the 
initial state l^);,^- Second, we turn off the couplings be- 
tween the charge qubits and the corresponding horizon- 
tal cavity modes and turn on those between the charge 
qubits and vertical cavity. We applies a coin flipping 
operation on the charge state. Third, we turn off the 
coupling between the charge qubits and vertical cavity 
and turn of those with horizontal ones and the free evo- 
lution of the interaction between the charge qubits and 
cavity modes is used to implement the conditional phase 
shift for one step of QW. 



IV. OPEN SYSTEM 

Coupling to additional uncontrollable degree of free- 
dom leads to energy relaxation and dephasing in the sys- 
tem. In the Born-Markov approximation, these effects 
can be characterized by a cavity photon leakage rate k 
and a pure dephasing rate 7 for each qubit of the coin 
state. The open system thus evolves according to 



^ = -t[JW] + Yl Ki*>[6j]P+ ifV[rt]p (41) 



with 



T>[L]p= -(2Lpfr - Dip- P DL). 



(42) 
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FIG. 2: (Color online.) The ln-ln plot of the phase spreads 
cr(4>) of the three-side coin QWs implemented in cavity QED 
with the DFT coin tosses and step size A8 — 0.8 as func- 
tions of the phase 4> f° r the initial coin states (|1) + i | — 1)) ® 
(|1) +i\ — l))/2 with the presence of decoherence. (a) Fixing 
the dephasing rate 7 = 0.06x, one can observe the QW-RW 
transition by increasing the cavity decay rate from k — to 
k = O.lx- (b) Fixing the decay rate k = O.Olx, the slope of 
ln-ln plot decreases slowly with 7/x increasing. 
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FIG. 3: (Color online.) With the step size A0 = 0.8, (a) 
AFDs and (b) ln-ln plot of the phase spreads a(4>) of the 
three-side coin QWs implemented in cavity QED with the 
DFT coin tosses as functions of the phase <f) for the initial 
coin states (|1) + i | — 1)) <8> (|1) + i \— 1)) /2 in the presence of 
decoherence. 



The master equation is used to compute p{t) from 
which the reduced state of the walker p w = Tr c p is ob- 
tained. As a signature of QWs the phase distribution can 
be obtained by performing full optical homodync tomog- 
raphy on the cavity to obtain the Wigner function, and 
hence the standard deviation thereby can be determined. 

With the realistic system parameters [3l|, |32[ 
(w a , W c , 5, e) /2tt = (7000, 5000, 100, 1000)MHz, simu- 
lated evolutions of the standard deviations of the phase 
distribution for the first several steps are presented in 
Fig. 2, which clearly reveal slope compatible with the 
characteristic quadratic decrease in phase spreading for 
increasing decoherence of the two-walker QW until the 
transition to the RW j36||. Here we use a DFT coin QW 
with initial coin state (|1) + i |-1)) <g> (|1) + i |-1)) /2 as 
an example with the decay rate of cavity k,\ = K2 increas- 
ing from to O.lx and the dephasing rate 71 = 72 = 
0.06x fixed. Thus the QW-RW transition is observed in 
Fig. 2(a). Furthermore, if the decay rate k\ = n^ = 0.01% 
is fixed, with the dephasing rate increasing the standard 
deviation of the phase distribution a as a function of the 
number of steps N decreases slowly. And with 7 increas- 
ing from to O.lx the slope of the In a- In TV plot decreases 
from 0.970 to 0.810 in Fig. 2(b). 

Except for the standard deviation of the probability 
distribution, one can determine the effects of the deco- 



herence on the quantum behavior of QWs via the AFD. 
As shown in Fig. 3(a), we have numerically calculated 
the AFD of the DFT coin QW with the initial coin 
state (|1) + i |-1)) <g> (|1) + i |-1» /2 in the two cases: 
(i) only with the decay of the cavity; (ii) only with the 
dephasing of the charge qubits. Using the parameters 
with which the standard deviations in the presence of 
the two sources of decoherence are quite close shown in 
Fig. 3(b), one can obverse that the AFDs of the two cases 
are quite different and the difference increases with the 
decoherence increasing. In Fig. 3(b), the \na-lnN plots 
are shown with different parameters. The result with 
«i = «2 = and 71 = 72 = 0.02x («i = «2 = 
and 71 = 72 = 0.2x) is quite similar with that with 
71 = 72 = and K\ — K2 = O.Olx (71 = 72 = and 
Ki = k-2 = O.lx)- Whereas, in Fig. 3(a) with the same 
choices of parameters, the AFDs are quite different and 
derive from each other. The AFD is more sensitive to the 
cavity decay rather than to the charge dephasing. These 
results imply the significant effect on the QW-RW transi- 
tion from the cavity decay k. Moreover, the cavity decay 
k is much more important than the charge dephasing 7 
with respect to the scaling of AFD with time t (propor- 
tional to the number of steps N) . The dephasing rate 7 
mainly leads to smearing the phase distribution and the 
phase distribution loses its symmetry. 



V. CONCLUSION 

In summary, we have introduced a protocol to imple- 
ment a QW with a multi-walker in phase space using 
cavity QED. The implementation of a QW with a multi- 
walker opens up the interesting possibility to introduce 
entanglement and more advanced walks. With different 
coin tosses and initial states the multi-walker QWs show 
different probability distributions which deviate strongly 
from the RWs and show faster spreadings and localiza- 
tion effect. By introducing decoherence, the transmission 
from QWs to the classical versions is observed. We pro- 
pose a physical realization to investigate the decoherence- 
induced irreversibility of QWs via the AFD. Our scheme 
provides an experimentally available way to monitor the 
detrimental influence on the QW by different decoherence 
sources. It is observed that the cavity decay has more 
detrimental effects on QWs rather than the dephasing of 
the charge qubits, which allows us to understand better 
the QW simulation in a realistic system. In conclusion, 
our theory establishes a pathway to realizing a many-step 
QW with multi- walker, and our techniques for observing 
the signature of QW via the AFD would be useful for 
general quantum information protocols. 
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